ORTHOGONALITY CONDITIONS AND ASYMPTOTIC STABILITY IN 
THE STEFAN PROBLEM WITH SURFACE TENSION 
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ABSTRACT. We prove nonlinear asymptotic stability of steady spheres in the 
two-phase Stefan problem with surface tension. Our method relies on the intro- 
duction of appropriate orthogonality conditions in conjunction with a high-order 
energy method. 



1. Introduction 

We are interested in the question of long-time nonlinear stability of steady state 
solutions to the two-phase Stefan problem with surface tension, one of the best 
known parabolic free boundary problems. It is a simple model of phase transitions 
in liquid-solid systems. 

Let Qcl" denote a C 1 -domain that contains a liquid and a solid separated by 
an interface T. As the melting or cooling take place the boundary moves and we 
are naturally led to a free boundary problem. Define the solid phase Q~(t) as 
a region encircled by T(t) and define the liquid phase Cl + (t) :=Cl\Q. . The un- 
knowns are the location of the interface {T(t);t >0} and the temperature function 
v : [0, T) X Q — > R. Let To be the initial position of the free boundary and vo : £2 — > R 
be the initial temperature. We denote the normal velocity of T by V and normal- 
ize it to be positive ifY is locally expanding 0. + {t). Furthermore, we denote the 
mean curvature of T by K. With these notations, (v,T) satisfies the following free 
boundary value problem: 

(1.1) d t v-Av = in Q.\F. 

(1.2) v=K on T, 

(1.3) V = [v n ]t on T, 

(1.4) v„=0 on d£l, 

(1.5) v(0,-)=vo; r(o)=r . 

Given v, we write v + and v~ for the restriction of v to Q + (t) and £l~(t), respec- 
tively. With this notation [v n ]t. stands for the jump of the normal derivatives across 
the interface T(t), namely [v n ]l :=v+— v~, where n stands for the unit normal on 
the hypersurface T(t) with respect to Q + (t). In (11.41) dQ. stands for the outer fixed 
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boundary of £2. Two basic identities related to the above problem are the "mass" 
conservation law: 



(1.6) d t \ / v(t,x)dx+\Q. (t 

and the energy dissipation law: 
rl 



= 0, 



- [ v 2 dx+\F(t)\\ + [ \Vv\ 2 dx = 0. 
2Jq. J Jq 



Here |£2~(f)| and \T(t)\ stand for the Lebesgue volume of Qr(t) and the surface 
area of F(t) respectively. Steady states of the above problem consist of static 
spheres and they form an (n+ 1) -dimensional family 

(1.7) ^:={l(/?,a)| aGM",/?GM + }, 
where for any a= (a 1 , . . . ,a n ) G W, R G M. + the pair 

E(i?,a):=((n-l)//?,S R (a)) 

is a time-independent solution of the problem (11.11 ) - (1 1 -4b . if Sr(a) Cfl. Let us 
parametrize the moving boundary T as a graph over a given steady state £(/?,a) G 
J^": we introduce the radius function r: [0,°o[xS" _1 — >M. and the parametrization 
0:[O,oo[xS"- 1 -»-r such that 

(1.8) 0(^)=a+r(f,§)§, ^GS"- 1 

and define the perturbation (u,f):=(y— (n— l)/R,r— 1). To (w,/) we associate a 
high-order energy norm given by 

N 

\\( u J)( t )\\ = Y,{W u \ I w^.-w^ao.fixfi) + 1 1 lw iV -^"iv 2 *+ 1 . 2 ([o,?] xS"- 1 ) } • 

k=0 

Our main result is the following: 

Theorem 1.1. Assume that C,r >0, where C, R := ^ — j^j^i fi.e. f/ze linear stability 
criterion holds). Then there exists e >0, smc/i f/ia? if initially 

||(«,/)(0)||<e, 

where (u,f)(t) := (v(f,-) — (n— 1)/R,r(t,-) — R) is the perturbation from the steady 
state £(/?,a), ?/ze?i ?/jere exijfcs a global-in-time unique solution (v,T)(t) to the prob- 
lem f li.il ) - rti.5D . Moreover, there exists a pair (i?,a) c/ore ?o (i?,a) smc/z ^/zaf (v,T) 
converges exponentially fast to £(^,a)GJ£\ More precisely, if we parametrize 
(v,T)(f) as a perturbation ofL(R,a), by setting T(t) = {a+r(f,§)§ , § G§" -1 } ant/ 
define (u,f)(t) = (v(t,-) — (n—l)/R,r(t,-)—R) then there exist constants C\ andC2 
such that 

||(«,/)WII<Cie- C2 ', t>0. 

Remark. By choosing iV large enough in the definition of the norm 1 1 • 1 1, we can 
make the solution as smooth as desired in a classical sense. We have not aimed 
for finding the minimal ,/V; rather, for the sake of clarity and conciseness of the 
estimates, we allow ourselves the flexibility of keeping ,/V sufficiently large. 
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Remark. Note that the time derivatives occurring in the expression ||(«,/)(0)|| 
are implicitly given by terms with only spatial derivatives, via the equations (11.11) 
and (Ob . 

1.1. Notation. On the unit sphere the Riemannian gradient with respect to 
the standard metric is denoted by V g and the Laplace-Beltrami operator by A g . For 
a given function h the k-th time derivative is interchangeably denoted by d t kh or 
hfk ■ When writing various norms of the functions defined on the unit sphere §"~ 1 , 
we drop the domain from the notation, e.g. | \f\ \ L 2 := | \f\ |z,2(gn-i). A ball (sphere) of 
radius R centered at a point a G R" will be denoted by fi«(a) (S#(a)). An important 
tool in our analysis are the spherical harmonics on the unit sphere. For a given 
function /jeL 2 (S' l_1 ), we introduce the spherical harmonics decomposition for h 

oo m(k) 

h = 22 H h k,iSk,i- 
k=0 i=l 

Here for each &gNq, stj, i=l,...,m(k) stand for the spherical harmonics of de- 
gree k. They are defined as restrictions of homogeneous polynomials of degree k 

in n variables onto the unit sphere. The set Uit=oUj=i f° rms an orthonor- 

mal basis on with respect to L 2 -product, thus justifying the above expansion. 
For more details, we refer the reader to li27l . The first n+l spherical harmonics 
*o,i)*i,i>---;*i,n will be denoted by so,s\,...,s n . For i=l,...,n we define the first 
momenta of h F\hi as well as the mean Poh by setting 

(1.9) mr.= [ hsidt;; F h = -^ T - [ hd$. 

We denote P2+h:=h— Pq^— Lfe=iPi^i- F° r a given function % we interchange- 
ably denote dj'W = % x i = e %fi and similarly for mixed partial derivatives. Einstein's 
summation convention is employed and we sum over repeated indices. A generic 
constant appearing in various estimates is denoted by C and it may change from 
line to line. The unit normal to a given surface T is denoted by «r and the index is 
dropped if the surface in question is clear from the context. 

1.2. Previous work. Stefan problem has been studied in a variety of mathematical 
literature over the past century (see for instance [34]). If (11.21 ) is replaced by the 
boundary condition v = on Y, the resulting problem is called the classical Stefan 
problem. It has been known that the classical Stefan problem admits unique global 
weak solutions in several dimensions ( lfI3Tl . lfl6l and [23]). The references to the 
regularity of weak solutions of the two-phase classical Stefan problem are, among 
others, Q, (H, 0, lfl2l . Local classical solutions are established in [22J and [26]. 

If the diffusion equation (11.11) is replaced by the elliptic equation Au = 0, then the 
resulting problem is called the Mullins-Sekerka problem (also known as the quasi- 
stationary Stefan problem or the Hele-Shaw problem with surface tension). Global 
existence for the two-phase quasi-stationary Stefan problem close to a sphere in 
two dimensions has been obtained in [8] and ifTTTl . and in arbitrary dimensions 
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in lfl4ll . Global stability for the one-phase quasi-stationary Stefan problem is es- 
tablished in |[T8l . Local-in-time solutions in parabolic Holder spaces in arbitrary 
dimensions are established in 0. 

Global stability of a single equilibrium for the Mullins-Sekerka problem follows 
also from the work of ALIKAKOS & FUSCO 0, as an immediate consequence 
of their study of Ostwald ripening. See the remarks at the end of Section Q] and 
Section 12. 1 [ respectively. 

In a recent preprint |[29l the authors analyze the thermodynamically consistent 
Stefan problem with surface tension. In the stable regime £#>0 their methods 
imply the stability of stationary spheres for the Stefan problem with surface tension 
using the maximal regularity techniques in L p -type parabolic Sobolev spaces. 

As to the Stefan problem with surface tension (also known as the Stefan problem 
with Gibbs-Thomson correction), global weak existence theory (without unique- 
ness) is developed in (31, ll24l and 13T1 . In ifTTl the authors consider the Stefan 
problem with small surface tension i.e. a <C 1 whereby (11.21) is replaced by v = OK. 
Local existence for the Stefan problem is studied in QUI . In lfl3l the authors prove 
a local existence and uniqueness result under a smallness assumption on the initial 
datum close to flat hypersurfaces. Linear stability and instability results for spheres 
are contained in ll28l . 

The first global-in-time nonlinear stability result for the flat steady hypersurfaces 
was given in [20l. Some of the references for the Stefan problem with surface 
tension and kinetic undercooling effects are iPTOl . ll30ll . 021 . ll33l . 

1.3. Motivation, methods and plan of the paper. To explain the linear stability 
criterion £r > and motivate our result, let us look at the linearization around a 
fixed steady state £(l,a) £ &\ 

(1.10) d t (u,f)=J?(u,f), 

where «£?(«,/) = (Au, — [d„u]'t), u = — {n — \)f — A g f on S n_1 . This linearization is 
easily obtained if we parametrize the moving boundary T as a graph over Si (a) : 
with the radius function r = 1 +/ and the parametrization as in (11.8I ). the pertur- 
bation («,/) takes the form (u,f) = (v—(n — l),R — l). It is readily checked that 
up to the first order Ko(j) = (n — 1) — (n — l)f — A g f and Vo$=f t . Associated to 
the linear problem dl.101 ) is the energy dissipation identity (cf. GT1 ): 

dt{ I u 2 + f {|V,/| 2 -(,i-l)/ 2 }} = -2 / |V«| 2 . 

Note that it is not even clear whether the above energy is positive definite, due to the 
presence of the negative definite term — Jg B -i (n—\)f 2 . Indeed, problem (11.101) may 
allow for a strictly positive eigenvalue under certain assumptions on the relative 
size of the domain £1 and the steady state L(R,a). More precisely, to any E(/?,a) € 
& we associate the stability parameter 

1 n-l 

Qr: ~\&\ \S r \R 2 ' 

where and \Sr\ stand for the Lebesgue volume of CI and the surface area of Sr 
respectively. It turns out that the sign of is critical to the stability properties of 
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L(R,a): if £#>0 the solution is linearly stable and if £#<0 it is linearly unstable 
(EH, [28]). The full nonlinear instability in the case £s<0, under the additional 
assumption that £2 is a perfect ball of a given radius R* is proved in [21 j. A closer 
look at the instability proof in |2~D shows that the proof itself is insensitive to the 
specific shape of £2 and works for general domains. On the other hand, in the 
stability regime (£#>0) the situation is more complicated. To explain this, note 
that the linearized problem (11.101) has (n+l) non-decaying solutions {oi},- = o,..., n , 
being exactly the (n+l) eigenvectors spanning the null-space of Jzf: 



where Sj,i=l,...,n, are the spherical harmonics defined in Subsection 11.11 The 
existence of a ( -s, i=0,...,n, encodes the (n+ 1) -dimensionality of the set & . We 
expect the perturbation of £(/?,a) to converge to a nearby asymptotic state L(R,sl). 
On the other hand, note that the first modes Pi/} = J s „-ifstd^Si (i = \,...,n) in 
the spherical-harmonics expansion of / are precisely cancelled by the expression 
Z(h) := /§,,-! { | V g h\ 2 — (n — l)h 2 } and they cannot be a-priori controlled by the en- 
ergy. This feature of the problem is another manifestation of the non-trivial null- 
space structure of the linearized operator «Sf . It causes a major difficulty in con- 
trolling the first momenta Pi/j-, i= 1, . . . ,n, thus introducing additional analytic dif- 
ficulties. 

In this paper, as stated in Theorem 11.11 we will prove the nonlinear asymptotic 
stability of the steady state L(R,a) € in the stability regime £r>0 for general 
domains £2. Our method has two basic ingredients. We start by introducing a set 
of geometrically motivated orthogonality conditions, that allow us to "mod out" 
an inherent degeneracy related to the existence of the non-trivial nullspace of the 
linearized operator. In analytical terms, we express (v,T) in a set of new, tubular 
coordinates. The second ingredient is the high-order energy method developed 
in EDI . ll2Tl as apart of the program to investigate stability and instability of steady 
states in phase transition phenomena. The goal is to prove an energy estimate of 
the form 



(1.11) g(t)+ / @(z)dT<£(0) + (S+C sup / 0(t)«/t, 



where $ and 3) are the energy and dissipation naturally associated to the problem 
(cf. (13.481) and (I3.49I )). If $ and 8 >0 are small, we can absorb the right-most term 
above into LHS, to obtain an a-priori estimate 



thus recovering the smallness assumption on £ if is small. Through an it- 
eration scheme and a continuity argument, this reasoning is made rigorous: proof 
of (11.111 ) is based on a series of energy estimates for the error terms. However, 
while in the previous works [20] and [21] the boundary conditions (11.21 ) and (11.31 ) 
honored the energy method fully, by introducing the tubular coordinates new error 
terms of (only) quadratic order emerge (Lemma I3T21 . The precise structure of the 
linearized curvature operator is exploited to get cancellation for such terms and to 



a :=(n- 1,-1), <7i 



(0,Si) i=l,...,n, 
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finally close the estimates (Theorem 15. lb . To obtain global existence for the solu- 
tion (v,r), we use a Poincare-inequality type bound $ <C3. It allows us to prove 
some decay of the solution on a bounded time interval. Together with a suitable 
smallness assumption on the initial data this decay can be bootstrapped to yield a 
global-in-time existence result. 

We wish to point out that our method can also be used to prove asymptotic sta- 
bility of steady spheres for the Mullins-Sekerka problem in arbitrary dimensions. 

Remark. After the publication of this work, the author learned that the sta- 
bility of Mullins-Sekerka steady sphere follows easily from the study of Ostwald 
ripening ||2] [3j where the problem of long-time behavior close to a steady state 
consisting of a finite union of steady spheres is analyzed. Moreover, the same 
orthogonality conditions as in this paper are introduced in (2J to deal with the 
non-trivial null-space - see the remark at the end of Section I27T1 

The plan of the paper is as follows. In Section [2] we heuristically motivate and 
then introduce the orthogonality conditions. We then derive the modulation equa- 
tions and reformulate the Stefan problem in the set of new, tubular coordinates. In 
Section |3l the Stefan problem is reformulated on a fixed domain and the high-order 
energy is introduced. In Section [4] we prove a Poincare type estimate (Lemma |4~TT ) 
and the positive-denniteness of the energy (Lemma |43T >. With these preparations, 
energy estimates are performed and local existence theorem is proved in Section [5] 
Finally, the proof of the main result (Theorem ll.il ) is presented in Section [6l 

2. Orthogonality conditions. Evolution problem 

To facilitate the analysis, without loss of generality, we assume a few simplifi- 
cations. We will perturb away from the steady state £(1,0) (assuming Si(0) c£2). 
Furthermore, note that the mass conservation law (11.61 ) necessarily determines 
the radius of the final asymptotic state. Recalling the parametrization (11.81 ) of 
F(t), dHU) takes the form d t M(y,T)(t) =0, where 

M(v,r)(t):= f v(t,x)dx+ f r -^^d$. 

We will henceforth assume that initially 

(2.12) M(v ,r )=M(I(l,0)). 

Condition (12.121 ) forces the asymptotic steady state to have the radius R aS ymp = 1 
due to the conservation of M{v,Y){t). This assumption constraints our stability 
analysis to the "manifold" of steady states ^ C & (defined in (11.71 )) consisting of 
the elements Z(a) :=Z(l,a) of the fixed radius R = l: 

S?:={£(a)| ae£l, 5i(a)cO}. 

Finally, we assume J^-i/o^f = 0, i = l,...,n where the initial surface Fq is 
parametrized by (l+/o(<^))<§, ^g§" _1 . Otherwise, translate the steady state to 
a nearby one, so that the condition is satisfied. 
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2.1. Heuristics. To motivate the analysis of the present work, in the following we 
provide some geometric heuristics for the choice of the above mentioned orthogo- 
nality conditions. Namely, let us think of <S as a submanifold of the subspace 
ofL 2 (n)xL 2 (S"" 1 ;n) of the functions of the form (u,&+R(£)£): 

JT:= {( Bj a+i-(£)§), aGM" t u:Q->R, riS"" 1 ^R} cL 2 (Q)xL 2 (§ n - l ;Q.). 

Thus, for aeR", E(a)(x,^) = (n— l,a+£) The non-decaying solutions 

{Oi}i=i,.,.,n of (11-1 Ob correspond to the infinitesimal changes in the center coor- 
dinate parameters at, i= 1 , . . . ,n. Motivated by this observation, we expect the so- 
lution to the full nonlinear problem to decompose into a component tangential to 
<S and the dissipating part that belongs to a plane transversal to £f in Jff. In fact, 
we shall demand that this plane is exactly the fiber ^- L , L 2 -orthogonal to in J4f. 
In other words, we choose the tubular coordinates (a(t),u(t,-),f(t,-)) in a neigh- 
borhood of the steady state manifold ( S: 

(2.13) (v(t),r(t))=-L a(t) + (u(t),f(t£)$), 
such that 

(2.14) (u(t,-),f(t,-)i(-))e&^, 

where i:S" _1 — >W 1 is the inclusion operator. Condition (12.141 ) is precisely the 
orthogonality condition. To obtain an analytic expression for (12.141) . we note that 
the tangent space T^p is spanned by the set {L i (x,E ! ) = d a , 1 L a(t) \ al=a , {t) } i= ^,_ n = 
{(0,e,)}, = i ... „. Thus, condition (12.141) implies that for any i= l,...,n: 

(2.15) 0=((u,fi),(0,e i )) L2 = f f(t,S)S- ei dS= [ f(t,Z) Si dS. 

Here e\ stands for the j-th unit vector in W. We have provided a geometric picture 
that renders the right orthogonality condition (12. 15b . thus drawing a parallel to the 
work of FRIESECKE & PEGO (HI on the stability of solitons in Fermi-Pasta-Ulam 
lattices. There, the authors introduce suitable orthogonality conditions exploiting 
the symplectic structure of the linearized problem. Their perturbed solution be- 
longs to the fiber symplectically orthogonal to the soliton state manifold. 

Remark. After the publication of this work, the author learned that the same 
orthogonality conditions as in (12.151) were first introduced and used by ALIKAKOS 
& FUSCO JH (see also Q) in the study of Ostwald ripening in the context of the 
related Mullins-Sekerka problem (see in particular eqns. (0.24) and (0.25) from flU 
and for the heuristic motivation see the remarks on top of the page 443 in |[2l). 
Their methods rely on optimal regularity techniques and, among other, imply the 
stability of a single equilibrium in the Mullins-Sekerka problem as pointed out in 
the introduction. Same idea is also used in a related work on the normalized mean 
curvature flow (T). 

2.2. Modulation equations. The evolution problem. A natural question raised 
by the form of the orthogonality condition (12.151 ) is to find the modulation equation 
(i.e. evolution equation) for the vector a(t). This is not straightforward because 
a(?) does not explicitly appear in (12- 15b . thus preventing us from differentiating in 



8 



MAHIR HADZIC 



time directly. To resolve this difficulty, for any i=l,...,n, let us define the mo- 
mentum test functions Pi(t,-) €C°°(0): pi(x,t)=x' —a'{t), where x=(x l ,...,x n ). 
Assume for the moment that (v,r)(t) is a classical solution of the Stefan prob- 
lem (11.11 ) - d 1.5b - We then multiply the equation (11.11 ) by the momentum test func- 
tions, use the integration by parts once and the boundary condition [d n v]t = Vy- As 
a result, for i= 1 ,...,n, we obtain the identity 

(2.16) / V T Pi= I u tPi + [ Vu-Vpi. 

Jv(t) Jo. Jn 

Note that the normal velocity expressed in local coordinates takes the form 

V r o0(^) = -^+a-H r of \ g \ = J r 2 + \V g r\ 2 
\g\ v 

and the mean curvature fc takes the form: 

Ko <j>(t,i;) = ^-]-V g -^ = (n-l)-(n-l)f-A g f+N(f), 
\8\ K \8\ 

where N(f) stands for the quadratic nonlinear remainder. The volume element 
dS(T) takes the form r" \g\d% in the local coordinates on the sphere, so we obtain 

/ V r pidS(r) = -f r"r^'+a-/ Pi n r = - f r n r t ^+ii-j V Pi 

Jr(t) Js"- [ Jr{t) Js"- [ Ja-(t) 

= -[ ^ l +^t)\or{t)\, 

where we used the Stokes theorem in the second equality and V/?, = e ; - in the last. 
Plugging this back into (12.161 ). we conclude 

(2.17) \Q-(t)\d'(t) = [ r n r t ^+ [ u tPi + [ u x ,. 

JS"-' Jq Jo. 

Thus, if (v,r) is the solution of the problem with the above choice of coordinate 
description of T, then a fortiori, the moving center components a'(t) satisfy the 
differential equation (12.171) . Written in terms of /, for any i = 1 , . . . ,n, the first term 
on RHS of dTTTT ) takes the form: 

f r»r tSi dt; = f (l+/)"/ f ^ = £ / ("V/to, 

where in the second equality we used the orthogonality condition (12.161) . implying 
in particular J s „- 1 f t S{ = 0, i = 1 , . . . n. The unknowns to be solved for, are the pertur- 
bation («,/) and the moving center a. Setting p = ... ,/?„), Stefan problem with 
surface tension (11.11 ) - (11.51) in the tubular coordinates (u(t, -),f(t, -),a(?)) takes the 



form 

(2.18) d t u-Au = in £2(f). 

(2.19) u = K-{n-\) on T(t), 

(2.20) [M„]lo0 = -^+ a -H r o0. on S" 
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(2.21) u n =0 on dQ,, 

(2.22) k (t) = -±-j^f (fif*f t Z + -* [ UtV +-±—[ Vu. 

I" { t )\k=iJs n - l \kJ |&2 {t)\Ja(t) | ft (01 Mo 

(2.23) / /£=0. 

(2.24) b(0,-)=mo; r(0)=r ; a(0) = 0. 

(2.25) M(l+ M() ,r )=M(£(l,0)). 

The idea is to exploit dissipative properties of the problem (12.181 ) - (12.211 ) to obtain 
a time-decay estimate for (u,f). Plugging that into (12.221 ). we hope to obtain a 
decay estimate for a. By bootstrapping this procedure, we will "drive" the solution 
to its asymptotic equilibrium. 

3. Fixing the domain and the energy 

We first describe the pull back of the problem ( 12.181 ) - ( 12.251) onto the fixed 
domain Q.\§ n ~ l and then define the high-order energies $ and Ql, deriving in 
particular the corresponding energy identities. From this point onwards, we will 
denote S e :=§ n ~~ . Define the following change of variables: 

®(t,x) = 7c(t,x)(x—a(t)), 

where % is a smooth scalar-valued function with the following property: 

f X _L) , |x-a(f)|-l<d, 

(3.26) n(t,x) = { r (^k ) d is small. 

{ 1, |x-a(?)| -1 >2J 

For any xGT, we may write X=a (0+K i*la(f)| ) !*-*(*) I =a 00+ r (£)£» by the def- 
inition of r. Hence 

^, . f , . x— a(t) , x— a(t) 1 x—a(t) x— a(t) 

|x-a(f)| l x ~ a (OI K4) l x ~ a (OI l x ~ a (OI 

Thus does the job for us: >Q\S^ . Note that this map is natural from the 

point of view of our geometric assumption that the evolving interfaces Y are, in 
fact, graphs over the unit sphere y. 

The inverse map x(x) to the above change of variables is given by 

x = a(t)+p(t,x)x. 

To convince ourselves that the function p(t,-):Cl\J7 > — »-R is well defined, we ob- 
serve that p has to satisfy the relation 

n{t,a{t) +p(t,x)x)p(x) -1=0. 

The existence of such a p can be established by the implicit function theo- 
rem, applied to the equation F(t,s,x) = 0, where F:IxRxR"->R is given by 
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F(t,s,v) = n(t,a+sv)s—l. Namely, F s (t,s,v) = (Vn-v)s+n >0 since 7T«1 and 
V7T = 0(V g /) is small. We define w:D.\y^R by setting 

(3.27) w(t,x);=u(t,a(t)+p(t,x)x). 

The heat operator d t —A on the domain II will transform into a more complicated 
operator in the new coordinates as stated in the following lemma: 

Lemma 3.1. The push forward of the operator d t — A with respect to the map 0: 
Q.— >y reads: 

(d,-A) # u = w,- aijW mi - biW? , 

where 

(3.28) a t j := K 2 Sij+2Kx j 7t x , +x l x j \VK\ 2 , hi :=2K x i +A7IX 1 - K t x l . 
Furthermore, 

2 

(3.29) k„]i = y(Hl)°4> 

The proof of this lemma is a straightforward calculation and it is presented in 
the appendix. In conclusion, on the fixed domain £L\y, equations (12.181 ) - (12.201 ) 
take the form 



(3.30) w t -a i jWij-b i w i = in £2, 

(3.31) w=Ko(f)-(n-\) ony, 

2 

(3.32) [ w l+= (_5£ + 4. nr o0) on^. 

\8\ \g\ 



The boundary condition (12.211 ) reads d n w = on d£l. The modulation equa- 
tion (12.221 ) is easily expressed in the fixed coordinates 
(3.33) 

w-wm&L © / ' / ^ + /n ( "" ,)o0|detz,e> " l+ /n v " oQ|detD0 " 1 )- 

where formulas (11.741 ) and (11.751 ) are used to express u t o® and Vmo@ in terms of 
the function w = uo@. The orthogonality condition (12.231 ) retains its form and the 
initial conditions take the form: 

(3.34) w(0,-)=w :=u Q o&; f(0,-)=f ; a(0)=0. 

3.1. The model problem. Let pt : Q. — > M+ be a smooth non-negative cut-off func- 
tion such that }l =0 close to 5? and \i = 1 close to dd and the origin Q. For any 
i= l,...,n let us define a differential operator 

D'=Ai^, + (l-At)^., 



By "close" we mean in an open neighborhood of prescribed positive thickness. 
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where d^i is the partial tangential differentiation operator defined in Cartesian co- 
ordinates through 

d?iu = d x iU— — Vw- — 
\ x \ \ x \ 

For given functions v,w G H l (£2) a simple integration by parts shows 

/ D l vw = - / vD'w- / vwv(jU)+ / vwn', 
Ja Ja Ja Jda 

where v(ju) := (d x i — d^)jj, = pjrVjU • A, thus giving us the integration-by-parts for- 
mula for the operator D'. For a given multi-index m = (mi,...,m n ) and a non- 
negative natural number s G No we define 

Df:=d t s d m . 

If we denote «Sf :=ajjd x i x j +bid x i the second order elliptic operator appearing in 
Lemma |3~T1 then we shall define the commutator 

P?,Jf\u:=E?S?u-SfEl?u. 

We formulate the following model problem, which is then used to derive the high- 
order energy identities. Corresponding to the equations (12.181 ) - (12.201 ) we analyze 
the equations: 

(3.35) %-aij%j-bi% = s& in £2, 

(3.36) < %o<t> = -(n-\)%-k g %+3§ ony, 

2 

(3.37) \d n ^}to(j) = ^-(-^ + a(t)-n r o(l)+^) ony. 

\g\ \g\ 

Lemma 3.2. The following energy identities hold: 
(1) 



(3.38) 

(2) 



X -d t J a P^ 2 +\dt J^{\y g X\ 2 -(n-l)X 2 }+ JjiV&AVV 

= I P+l Q+ I <%%C, 
Ja Jy Jda 



p% 2 + l-d, ( p\DW\ 2 + d, f pV'% t AV'% + f pa ii D k %D k '%r j 
a 2 Ja Ja J a 

(3-39) +^{|V^| 2 -(n-l)| Zf | 2 } + l5 f ^{|V^| 2 -(„-l)|V^| 2 } 

= / S+ f T+ f [-(%+AW)W n +VW-VW n ), 
Ja Jy Jda 

where 

(3.40) 

P:=-{a ij p) xJ '^' i ^ + +bi%)p^\ Q:=^ah)-n r oQty 

r r 
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S := (jS+bM) {% ~ pDiDiW) + ^p,\DW\ 2 + {paij) t %% 

(3 " 41) -(aijp)j%% + (a ijP )j%D k D k ty -D k (a ijP )%D k Wj 

-pDjty %v{\i)- a u p %D k <%j\{\l) ; 

(3.42) 

T:=-Zt&t+— ^KW+^a(t)-nro<b%-^a(t)-nro<i>K !i W. 

r r \g\ \g\ 

Proof. We multiply ( 13.35D by p% and integrate over Q.. Integrating by parts, we 
obtain 

\d t [ pfy 2 + f pV^'AWty- f aij[%}n } pfy 
2 Ja Ja Jy 

= -[ p t <% 2 - [ (a ij p) xi W i ®r+ f {^+bi%)p^+ f 
2Ja Ja Ja J da 

Note that 
(3.43) 

- / a ij [%]n i p ( % = - I a ij n i n j r[^f n }'^ 

^ 2 2 

= "/^^ r R ( "t[ +a( °- HrO0+ ^ ( " ( ' 1 " 1)z " A ^ + ^ ) 

2 Jy Jy r Jy r 

This proves the first identity. As to the second identity, we multiply (13.351 ) by 

p 0W t —DiDffl^, integrate over Q. and use integration by parts: 

(3.44) 

I P% 2 + ~d t f p\DW\ 2 + d t f pV% t AVW + f p aii D k %D k q/; 
Jo. 2 Ja Ja Ja 

- f a ij n i n j r[%]%+ f a ij n i n j r[^ n ]tA g ^ = \ (srf +b i %)[%-pD i D i '%) 

+ ^ [ p t \DW\ 2 - f pD i W%v(ii)+ f {p aii ) t %Wj- f { aij p)j%% 
Ua Ja J a Ja 

+ / ( ai jp)j%D k D k W- [ D k {a ijP )%D k %- [ a ijP %D k %v{ll) 
J a J a Ja 

+ / [-(^+A<2r)<2$,+V<2r-V<2d. 
J da 

Let us extract the energy contribution from the integrals over 5? in the above iden- 
tity. 

- J a^nh^r^t = - M^r^ ( - ^ + a(t) -H r o<j> +ff ) ( - („- l) Xt - A g Xt +& t ) 



ly r A \g\ \g\ 
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Similarly, 

I 1 2 2 

[ tHjrtn'riVAtAaW = f H~r-- (- yr + a{t) ■n r o<j>+^) (- (n- l)A g x-A 2 g X+*g&) 

= \dt [ {\V 2 g X\ 2 -(n-l)\V gX \ 2 }+ [ r fMt)-n Y o<t>A g W+ f X A g @-—^ g W. 

From the two previous identities we obtain the expression (13.421 ) for the eiTor term 
T. □ 



3.2. The energy. We introduce a parameter y> 1 to be fixed later. Motivated by 
Lemma l3~2l for given 'i'j we define the auxiliary energy 



2Ja 2Ja J a 

Jj\V g X\ 2 -(n-l)X 2 } + yj\Vlx\ 2 -(n-l)\V g X\ 2 } 



(3.45) 



and the auxiliary dissipation 
(3.46) 

9(W,Z)-=y[ f p% 2 + f pa ij D k %D k W j + f {\V gXt \ 2 -(n-\)\x t \ 2 }. 

Jo. Jo. Jo. Jy 



Summing the two previous expression and using Lemma [3721 we obtain 
(3.47) 

£(W,Z){t)+ f 9(W,z)(s)ds=*(W,x)(P)+ f I {yP+S}+ f f {jQ+T} 

+ 11 \ a K(r^-^ t -A^)+y^-y^n\- 

Jo J dQ. 

We define the total energy and dissipation by setting 
(3.48) <?(*,/)= £ <?(Dfw,DTf), 



\m\+2s<2N 



(3.49) @(w,f)= £ @(D™w,D?f). 

\m\+2s<2N 

We shall often write S{t) :=S(w(t)J{t)) and 9{t) :=9){w{t)J{t)). 



4. A-PRIORI ESTIMATES 

Lemma 4.1 (Poincare type estimate). There exists a positive constant jS such that 
on the time interval of existence of a solution to the Stefan problem, the following 
estimate holds: 

\\f\\L2 + \\w\\ L 2 m <P\\Vw\\ L 2 {a) . 

Proof. To prove the lemma, it is much more instructive to work with the solution 
// on the moving domain Gl(t) and later convert it into a result for w. Recall the 
expansion: uo(f) = — (n—l)f—A g f+N(f). Fix a spherical harmonic s^i of degree 
k, where k > 2 and 1 < i<N(k) (recall the notation introduced in the introduction). 
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Multiply both sides of the above relation by s^j and integrate over 5? . Observing 
that J> (- (n - 1 )f-A g f)s u = { n -\)(k 2 -\ )f Li , we get 



(4.50) 
Note that 



fk,i 



1 



(n-l)(k 2 -l)Jy 



1 



(n-l)(k 2 -l)J,y 



N(f)s k ,i. 



y 



uoQsu 



1 



uod) 

y v \y\ Jy 



UO 



- C H M0 ^~ Tl^ J y uo $\\L2 <C\\ Vu \\lHQ±)-, 



\y\Jy 

where we used the Sobolev inequality in the last estimate above. The last inequal- 
ity, together with ( 14.501 ), immediately implies 



(4.51) \fkj\<cV&+C f N(f)s u d% 

Jy 

In order to estimate Jyf, we observe that 



k>2 



^- 2 \g\ = -{n-\)f-A g f+q{f), 



where q stands for the nonlinear remainder with a leading order quadratic term. 
Integrating the above equation over 5?, we find 



(4.52) 



(n-l) /+/ q(f), 



y 



y 



Multiplying the conservation law Jq, u + lyf +HL2H) lyf / n = () bv pj 

and ( 14.521 ) by and subtracting the two equations, we obtain 

(4.53) 



1 



1 



Cl\Ja \T(t)\Jr{t) ! n\a 



1 



k=2 



y 

1 n-l 



f k + 



1 



<?(/)• 



Note that C-Cp-f^)^^-^ and hence ^-^-^ 



\<C 



which for small enough implies | ^ — if^l > C/2 >0. Hence, upon di- 

viding (14.531) by Ki := ^ - jf^jy), we conclude 
(4.54) 



1 



^ V |^| 7^ \T{t)\Jr{t] 

From the mean value theorem, we deduce 



1 



1 



nKi\a\£ 2 



<C\\Vu\ 



y 



f k + 



1 



Ki m Jy 



1(f)- 



\Q\ Jn \r(t. 
Thus, from (14.541 ) and the previous inequality, we have 



L 2 (£2±)- 



(4.55) 



|^/|<C||V M || L2(fl±) +C||/||2 2 
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Summing (14.551 ) and (14.511) and keeping in mind that Pi/ = 0, we obtain 

. OO N(k) OO f f 

l/^/l + I I l/Ml<C||V M || i2(£J) +C||/||2 2 +c££{ J y N(f)s u 
<C||V M || L 2 ( ^+C||/||2 2 +C||^(/)|| i2 <C||V M || L 2 (fi) +C||/||2 2 . 

Smallness of \\f\\ii then implies 

oo N(k) 

L IA'l< c ll v "llL 2 (nw- 

k=0 i=l 

Changing variables (x— >&(x), i.e. u— >w) and using the smallness of ||V g /|| £ 2 we 
get the bound ||Vm|| L 2(^±-) < C| | Vw| |z, 2 (si)- Combined with previous estimate we 
conclude 

(4-56) ll/lb<C||Vw|| i2(£1(f)) 

Finally, from the previous inequality and the conservation law (11.6I ). we immedi- 
ately deduce 

I u 



<C||Vii|| L 2( Q ( f ))<C||Vw|| L 2( Q(f )). 
Changing variables and using the smallness of | \V g f\ | £2 , we obtain 



(4.57) 



w 

Q, 



<C\\Vw\ 



L 2 {Q 



and this finishes the proof of the lemma. □ 
In the new coordinates, the mass conservation law dl.61 ) takes a different form. 
The following lemma expresses this conservation in a way that will be useful for 
proving the positive definitiveness of the energy expressions $ and 

Lemma 4.2 (Mass conservation law). The following identity holds: 
(4-58) / pw+ [ /=- f £ (l)-+ I wg(Vp,f), 

where g is a bounded smooth function with g(0,0) =0. 

Proof. Integrating the equation ( 12.181 ) on a moving domain £2\T(?), using the 
Stokes formula and the boundary condition (12.201) . we obtain: 

3,1 „Wfi±QU. 

Ja(t) Jy n 
The assumed initial condition (12.251 ) implies 

J off) k = x W 

To express the integral over £l(t) above as an integral over Q\y , it remains 
to understand the Jacobian |detD0 _1 |, where 0:£2(f)— >Q\y is the change 
of variables map defined in the line before (13.261 ) and 1 :x — > a(t) + p (t ,x)x. 
Thus, {D<d~ { )ij = p8ij+Px>x ] where 5,-y denotes the Kronecker delta. This im- 
plies detD© 1 =p n +qi(Vp) = p + (p — \)q2(p)+qi(Vp), for some polynomials 
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qi and q 2 such that g,(0) =0, i= 1,2. From here, we easily infer the lemma. Since 
p — 1 =/ close to y and p is smooth and equal to close to the boundary d£l, the 
claim of the lemma follows. □ 

Lemma 4.3 (Positivity of the energy). Under the smallness assumption on 
\\D™ f\\ L i + \\D™w\\ L 2, for \m\+2s<N, the energy quantities <§ and 3) are posi- 
tive definite. 

Proof. For a given function co :^— >R, let us abbreviate Z(co) = fy\V g co\ 2 — (n — 
I) CO 2 . If %=T)™f, with |;ll|>l, then, by Wirtinger's inequality, we immediately 

see 

Z(x)>C\\F 2+ x\\ 2 L2 

and analogous inequalities hold for Z(% t ) and Z(V^). If however, % is of the form 

X = d t sf, then we must exploit the conservation law ( I4.58I ). Note that 

(4.59) 

/ a ^ +z ^ ) =4 ( /n pw ' )2+ / a p(w ''"^X p ^ )2 "w ( // , ' )2+z(p/ '' 

where \Q\ p = f n pdx. From Lemma 14. 21 we obtain 



/ pw,*+ / fs=G(w,f), 
Ja Jy 



where G(w,f) = f a { -d t * (pw) +pw t *}-f y ZU ©/rV«+/^(>vg(Vp,/)). 
Thus, 

= ^-W )( /n PB ''' )2 - G(B '' /)(G(W>/> - 2 /n PB '' ,) - 
Note however 

1 n—\ y 1 1 y , , , r \ 

w P ~w\ =c+ wrw P =c+ 

where h is a smooth bounded function, /z(0) = 0. Going back to (14.591) . we conclude 

f pw 2 s+Z(f t ,) = C( f pw t s) 2 + f p( Wt ,-JL f pw t s) 2 +Z{Ff tS )+R( Wt s,f tS ) 
Ja Ja Ja \Ll\pJa 

where the remainder R(w t s,f t s) is a cubic nonlinearity given by 

R(w t s,f t *):=-G(w,f)(G(w,f)-2j^pw tS )+h(f)(j^pw tS ) 2 . 

Under a smallness assumption on w and / it is easy to see that .R(wv,/r s ) < 
Mll w Hlz2(n)~HI^Hli2), for a small constant A and thus there exists a positive 
constant M > 1 such that 

^ (I | ||z (£J) + 1 1 1| 2 ) < +^(/^) <m(| | |i 2(£2) + 1 1/,, I . 

This concludes the proof of the lemma. □ 
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The following lemma states that the norms defined by $ and *2l are equivalent 
to the corresponding parabolic Sobolev norms. The main new ingredient is that 
all the derivatives of w are bounded. In particular, derivatives of w in the normal 
direction are bounded by the energy quantities S and $l, although they are not 
a-priori contained in their definitions (13.481) and (13.491) . 

Lemma 4.4. The following norms are equivalent: 

N 

{ W J) (0 ~ { 1 1 W 1 1 W N - k -""W 2k ' 2 ( [0.f] x £2) + 1 1 ^gf\ I H"»-'- W 2k + 1 2 ( [0,1] x S»- 1 ) } • 
£=0 

N 

(0 ~ { 1 1 V>V I lw w -*. 2 W 2 *+ 1 . 2 ([0,f]x«) + 1 1 W ' 1 1 W*-tV,2([0,t]xfl) 
fe=0 

+ 11/* I I ^W-*,2 W 2i+l,2([0 jf ] xfl) } ! 

where A ~ B means that there is a constant c > 1 smc/i f/zaf - < B < cA. 

Proof. We only sketch the proof as it is analogous to the proof of Lemma 3.5 
in EDI . What we need to prove is that for any triple of indices {}JL,r,s) satisfying 
\{J.\+r+2s<2N we have 

\\I%d&M\li<n)<C*{w,f). 

The claim is evidently true for r = 0, 1 following from the definition of S. We then 
proceed by induction on the number of normal derivatives r. Inductive step uses 
the basic equation (13.30I ). It allows to express w nn as the sum of the multiples of w t , 
D'^w and D'w n in the region of Q. close to the unit sphere 5?, thus allowing to com- 
plete the proof in the case r = 2. By successively differentiating the equation (13.301 ) 
and expressing the term with the highest number of normal derivatives in terms 
of terms with either less normal derivatives or purely tangential derivatives, we 
inductively complete the argument. □ 



5. Local existence 

Theorem 5.1. There exist small positive constants E\, m\ such that for any Eq < E\, 
mQ<m\ and positive constants T*, C*, such that if 

^(w ,/o)<£o; |a(0)|<iiio, 

then there exists a unique solution to the Stefan problem A2.18i - (12. 251 ) on the time 
interval [0, T*] such that for any < s < t < T*: 

(5.60) sup <f(r)+ / @(%)dT;<£{s) + {\+C sup y/S{x)) [ 9{x)dr; 

(5.61) sup \a(i)\<2mQ. 

0<r<r* 

Moreover, S is continuous on [0,T* [ and sup ?<T< , 
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Proof. We first sketch the proof of the a-priori bound (15.601 ) assuming that the 
solution already exists. We explain in detail how to bound the hardest quadratic-in- 
order terms. For all the other (cubic) terms we use rather standard energy estimates, 
that can be found in the corresponding local existence proofs in [20] (or [21]). 
To prove the estimate (15.601) . we must bound the right-hand side of the energy 
identity (13.471 ) where we plug in =Dfw and % =Dff for all \m\ +2s< 2N. All 
the terms in the definitions of P, Q, S, and T (cf. (13.401) - (13.421) ) are trilinear except 
for the underlined terms, which are only quadratic in the order of nonlinearity and 
it is a-priori not clear how to bound them. 

There are two types of the quadratic error terms: the first type arises for technical 
reasons due to the introduction of the cut-off function /I while fixing the domain in 
Section [3] These are the two underlined terms in the expression (13.411 ) for S. For 
% =D™w and %=D™f, \m\ +2s<2N, by the Young's inequality, we have 



(5.62) 



/ o 7 n p Di **,v(,)|<I|'/ n P^ + ||v(,)||i„, a) // n P|D*r 
4 Jo J a Jo Jo. 



Similarly we estimate the integral of the second underlined term in the expres- 
sion (13.411 ): 
(5.63) 

/ / a ijP %D k %v{pL) <- [ [ a ijP D k %D u %+C* f f pVW'AVW , 
Jo Jo, 4 Jo Jo. Jo J a 

where we possibly enlarge the constant C*. The size of the constant C* dictates 
the choice of y in the definition of the energy (cf. Section 13^21 : set y=3C*. We 
can thus absorb the right-most terms in the estimates (15.621 ) and (15.631 ) into the 
/-dependent term in the definition (13.461) of Qi. 

The second type of quadratic error terms is completely intrinsic and arises due 
to the presence of the moving center coordinate a(?) in the parametrization of the 
moving surface T{t). They are the three terms, underlined in the expressions (13.401) 
and (13.421 ) for Q and T respectively. Note that, if \m\ > 1, then a(t) =D™a(t) =0, 
since a depends only on t. Thus, we are only concerned with estimating expres- 
sions containing a(t) of the form a(t) = d t k+l a(t) = a^ +1 ' (t), where 0<k<N — 1. 
The underlined expression in the formula for Q thus takes the form 



f ( ^ k+i \t)-n r o<pdl i w. 
Jo Jy \g\ 



The difficulty is immediately clear: since the expression ^|Hro0 is a quantity of 
order 1 for small /, the whole integral is only of a quadratic order and it is hence 
unclear how to bound it by v© JoS*- In the local coordinates on the sphere: 

n r o<f) = —%--±-, 
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where we keep in mind that £ is the unit normal of the unit sphere 5?. From here, 

" 2 - l g\\ rV g f 



Using this identity, 



o Jy 



a { - k+l >(t)-n r o(f)d t k w-- 



Jk+l) 



o Jy 



d*(-(n-l)f-A g f+N(f)) 



The key observation is t, GNull(A^ + (?i — where J stands for the identity 
operator. Hence 

^a^ +1 )(0-^ f *(A, + ( ? 2-l)^)/ = a(' :+1 )(0-^(A, + (n-l)^)^^ = 0. 
Using the previous two identites, we eliminate the purely quadratic contribution: 



o Jy 



a^ +1 )(0-[^ + - 



,(*+!) 



o Jy 



rV f 

—r^]dA-(n-l)f-A g f+N(f)) 
rV f 

-^f](-(n-l)f tk -A g f tk ) 



o Jy 



i {k+l Ht)-[^- r ^f]dAN(f)) 



Recall that N(f) is a quadratic nonlinearity and therefore the two integrands on 
the right-hand side above have a manifestly trilinear structure and are thus easy 
to estimate. Namely, applying the differential operator d t k to (12.221 ). it is easy 
to deduce the bound |a^ +1 )(?)| <C^/~3l if we knew that |a(T)| is bounded. To 
estimate |a(r)| we use (12.221) . Recalling p = x— a, and assuming smallness of S, 
we obtain 



|a(r)|<|a(0)| + r sup |a(T)|<|a(0)|+C7V<f+CT sup |a(T)|||K f || L 2 (Q) . 

o<-r<r o<r<r 

Therefore 

sup |a(T)|(l-C7V^)<|a(0)|+Crv^, 

o<t<r 

thus implying sup 0<T<r |a(T)| <2|a(0)| <2/«o f° r appropriately small S" and some 
finite T . Integration by parts, Lemma |4~T1 and standard energy estimates then imply 

h Jy \g\ Jo 
The same idea as in the proof of (15.641 ) works for the remaining two underlined 
integrals appearing in JgJyT (cf. (13.42I )). so we finally conclude (for a=a^ k+1 ' (t) 
and =w t k+i) 



(5.64) 



(5.65) 



Jo Jy \g\ Jo Jy \g\ Jo 



The construction of the solution (w,f) follows identically the construction scheme 
for the local solution from ll20l . We briefly summarize the main steps: we set-up an 
iteration scheme, which generates a sequence of iterates {(w m ,/ m )} me N solving a 
sequence of linear parabolic problems. As in [20], such an iteration is well defined, 
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but it breaks the natural energy setting due to the lack of exact cancellations in the 
presence of cross-terms. We design the elliptic regularization 

\g\ r \g\ 

to overcome this difficulty. For a fixed e, we use it to prove that {(w m ,/ m )} m6 N 
is a Cauchy sequence in the energy space. Upon passing to the limit m— >oo we 
obtain a solution existing on an a-priori £ -dependent time interval [0,r e ]. But, in 
the limit m — } °° the dangerous cross-terms vanish. Since the elliptic regularization 
honors the energy structure, we obtain an energy bound analogous to (15.601 ) with 
e-independent coefficients. By continuity, the solution exists on an £-indpendent 
time interval [0,T]. Finally, we pass to the limit as e— >0 to obtain the solution of 
the original problem. □ 

6. Global existence and asymptotic stability 

The proof of the main result Theorem 11.11 is an immediate consequence of the 
following theorem and Lemma |4~4l about the equivalence of norms. 

Theorem 6.1. There exist small positive constants E and m, such that if 

£(w Jo)<E; |a(0)|<m, 

then there exists a unique global-in-time solution to the Stefan problem < \2.18[ - 
< \2.25i converging asymptotically to some steady state solution £(a), where 

a(f)— >2L ast—>°°. 

Moreover, there are constants C\,C2 >0 such that the following exponential decay 
estimate holds: 

£{w{t)J(t))<c x e- c «. 

Before we prove the theorem, we will prove an important auxiliary estimate, 
necessary for the proof of Theorem l6.ll 

Lemma 6.2. Let [0,t] be an interval of existence of solution to the Stefan prob- 
lem < \2. 18i - (12. 251 ) for which the estimates H5.60\) and H5.61\) hold. Then there exist 
constants a, 8 Gl + such that 

„, , 40<f (0) at . 
g(t)< - t y ' e~-; t€[0,t[ 

and 



(6.66) m < c J^Ee-9 ; te[0,t{, 

\/t 

where j8 >0 is given by Lemma \4A\ and cp§ :=2a//3c). 

Proof. As in [25 ], p. 135, we define V{s) := f^(x) dx. From (I5T601) and 
Lemma |4~T1 we conclude that there exists a constant a >0 such that 

(6.67) sup £(%) + a [ £{%)dT<£{s). 

s<r<t Js 
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Thus aV{s)<£{s). On the other hand, V'(s) = -£(s)< -aV(s), which, in turn, 
implies: V(s) < V(0)e~ as . Integrate (16.671 ) over the interval [f/2,f] with respect to 
s, to obtain: 

^sup^(T)<V(i)<V(0) e -f, 



implying 



/ x 2V(0) m 
(t)< — —e~-. 



Note that for Eq and mq small enough (g (0) < 3^2) in Theorem O by continu- 
ity and the inequality (15.601 ). we obtain the bound Csup 0<s<t £{s) < \. Plug- 
ging back into (15.601 ) and absorbing the right-most term into LHS, we obtain 
J2#(t)«Zt<2<£(0) for all fe[0,t[. Hence, by LemmagT] V(0) = / V (t)Jt< 
PfQ®(z)d'V<2pg' (0) and we obtain 

(6.68) /( <)s a^-f, 

thus proving the first claim of the lemma. From (12.221 ) and Jensen's inequality, 
we conclude that for some constant 8 >0 on the time interval [0,t] the following 
inequality holds: 

8 n 8 8 

|a| < jLll/lbll/fll^ + jIKIIia^ + jll^Iha) 
^ k=i * * 

* k=\ 



By (16.681) we obtain 



for any f € [0,t] and this finishes the proof of the lemma. □ 
Proof of Theorem 16.11 Let m = ^ (mi is given by Theorem 15. II ) and E be such 

_ ay 

that cp§c a y/E< whereby Ca :=max v >7-* ^ ^ . Define 

(l-e 3~) 

(6.69) ^:=sup{ sup |a(r)| <2m& £20) holds for any s,t G [0,t]). 

Assume ^<oo. With T* >0 given by Theorem 15. II choose a unique LgN such 
that 

^e[Lr*,(L+i)r*[. 
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Integrating over [kT* ,(k+l)T*] with l<k<L— l,we obtain 
|a((fc+l)r*)|| < |a(itr)|+ [ {k+l)T \a(x)\dx 

JkT* 

(k+l)T* 



< \ a (kT*)\+c pS y^) 



e 4 



kT 



< |a(*r*)| + 



'kT* 



e 4 dx 



kT" 



<\ a (kn\+ cps r^ T*e-^ 



'kT* 



< \a(kT*)\+ Cl}S ^I{0)Vf^e- 
Summing over k=l,...,L— 1, we obtain 



akT* 



L-1 



i(LT*)\ < \*(T*)\+cpsy/£(ti)y/T*£e 



akT* , . _ , 

— = |a(r*)| + 



cpsVW)VT*~ _«£ 

— : \P. 4 



k=\ 



l-e 



< 2m ■ 



CrsVT* _gr^ . — , — 3 

—e 4 \/E<2m+cp S c a vE< - 



l — e 4 



where the last inequality follows from the choice of m and E above. If we denote 



-LT* 



,9 -IT* 



we obtain 



3mi 



\a(T l )\ < |a(Lr*)|+ / | a (T)|jT<— i+c^T^O)^- 
J IT* 4 



aLT"- 
4 



<— +cp S c a VE<— + — 



-mi. 



Since ^(r 1 ) <£ <£i (with E possibly smaller, cf. (I6T68T )) and ja^ 1 )] <mi, we 
can extend the solution (w,/) starting at ? = T 1 to the time interval {0,T l +T* [, by 
Theorem 15. II (local existence theorem). Same theorem guarantees the estimate 

ll a llz.°°[o,r 1 +r*] <2mi 

as well as the validity of the energy estimate (15.601 ) on the time-interval [0,T + 
T*]. Since we assumed < °°, we obtain T 1 + T* > 2F and this, together with the 
continuity of £(w,f) and a contradicts the maximality of . It is now evident that 
the decay estimates (16.661 ) and (16.681) from Lemma [6721 hold for all t 6 [0,°°[. This 
implies the decay claim of Theorem 16.11 as well as the existence of an asymptotic 
center aGll such that a(f) — >a, since a(t)— >0 as t— >oo. □ 



Appendix A. Proof of Lemma |3~T1 
Note that (13.261 ) implies 
( 1 .70) x x i = nej + % x i (x — a) . 
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From n(x) = l/p(x), we obtain K x i = — (yp-x x i)/p 2 , which in turn, combined 
with (11.701 ). after an elementary calculation implies 

where 

*¥=p+Vpx. 

After further differentiating (11.711 ) with respect to x l and using the relation (11.701 ). 
we arrive at 

-Ap p x , xJ p x ^{\-p 2 ) -p^¥|Vp| 2 +2|Vp| 2 Qy+p) 

p 3x i> p 5x i> 2 p 3x F 3 

Similarly it is not hard to see that n, = (— p, +p Vp -a) / p*F. In order to evaluate the 
Laplacian in new coordinates, by ( 13.271 ) we first write 

(1.72) u(t,x)=w(t,7i(t,x)x). 
Applying A x to the left-hand side, we obtain 

(1.73) Au(t,x) = n 2 A x w+ (2nx J n x i +x i x i \Vn\ 2 )w x i X j + [l% x i +Anx')w x i 
Applying d t to the left hand side of (11.72I ). we obtain 

(1.74) Uf(t,x) = w t + K t V x w-x. 
Since (d t — A)u = 0, we conclude from (11.731 ) and dl.741 ) 

(d t - A. f )w = dijWjig +bi\v x , , 

where a,j and bj are given by (13-281 )- Furthermore, using d 1.7 II ), it is easy to see 
that 

-_L* ^ . *'* ; 'I v pI 2 

We now turn to the proof of (13.291 )- Note that 



u x i = V x w ■ x x i = V x w ■(7ie i + 7i x i(x-a(t)))= V x w • ( ^ - px) . 
From here, we infer the formula 

(1.75) v^H-Vwf;)-^^^). 

From the above formula, we obtain 

(1.76) [vu}± °m=-^ « ) - ) g • [ y H - (g ) 

Note that 

5 • [v w ]i(0=£ -([vrf+M = Mi> 

since [V g w]i=0 (due to the fact that u + \y = u~\y) and ^ -ny = \ny\ 2 = 1 (ny 
stands for the unit normal on S^). Observe further that p(£) =r(£). Furthermore, 
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since Vp(£)-£ =V g p-£ =0, we also have x V(t;) = p(t;) = r(t;). Form these obser- 
vations and from (11.761) we obtain the formula 



[V//] c <;) 1 [„■„] W - [w n ]t = - r [w n }tny - -f [w n ] _ . 



It is straightforward to see that n' r o = — Vg ^J s ^ . Hence 

[V»]± -n r ) o0 = (l[w„]±^- ^f^' K]l) • [j-^- Vgf ^' 

= ^W«]t+ l ^[w n ]±= l 4[w n ]t 
\g\ \g\r l r l 

and this proves ( I3.29I ). □ 

The author wishes to thank Yan Guo for many stimulating discussions. This 
work has been partly supported by the NSF grant CMG 0530862. 
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